Local zeta regularization and the Casimir 

effect 



In this paper, whose ahns are mainly pedagogical, we illustrate how to use 
the local zeta regularization to compute the stress-energy tensor of the Casimir 
effect. Our attention is devoted to the case of a neutral, massless scalar field 
in flat space-time, on a space domain with suitable (e.g., Dirichlet) boundary 
conditions. After a simple outline of the local zeta method, we exemplify it in 
the typical case of a field between two parallel plates, or outside them. The 
results are shown to agree with the ones obtained by more popular methods, 
such as point splitting regularization. In comparison with these alternative 
methods, local zeta regularization has the advantage to give directly finite 
results via analitic continuation, with no need to remove or subtract divergent 
quantities. 
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1 Introduction 



Zeta regularization is a method to give meaning to the divergent series appearing 
frequently in quantum field theory, reinterpreting them as analytic continuations. 
For example, the divergent series 

oo 

E^' (1.1) 

£=1 

is interpreted in this approach as the analytic continuation at s = —3 of the regu- 
larized series ({s) := Y2t=i ^^at converges for 3fJs > 1 and defines the familiar 
Riemann zeta function; in this sense, the sum (11. ip "equals" C(~3) = 1/120. Tricks 
of this kind have been used for a long time in quantum field theory: for example, 
the above analysis of the series (II. ip appears in one of the most popular derivations 
of the total Casimir energy for a scalar or electromagnetic field between two parallel 
plates (see PQ HI E], or the issue "Casimir effect" in Wikipedia). 
The computation of local quantities, such as (the vacuum expectation value of) the 
stress-energy tensor, can be done as well via a generalization of the above method; 
this procedure, called the local zeta regularization, is a bit less popular than its 
analogue for global quantities, such as the total energy. 

The method of local zeta regularization arose from some ideas of Hawking [6] and 
Wald [H]; these were systematically developed and applied to the stress-energy 
tensor by Moretti in a long series of papers, among which we quote [9l [TO] . 
These authors typically work on curved space-times, in a Euclidean framework (i.e., 
with a space-time metric of signature (-I-, -|-, -|-, -|-)) (indeed, [6] |9] also mention a 
flate case, namely, a four dimensional Euclidean torus); in this setting, the local zeta 
regularization is applied to divergent sums arising from path integrals. 
Since the local zeta method is not so popular, in our opinion it is not useless, at 
least for pedagogical reasons, to illustrate it in a much simpler framework; this is 
the aim of the present work. 

In this paper we consider a (neutral, massless) scalar field in Minkowski space- 
time (with the usual metric of signature (— , +, +, +)), as viewed in a given inertial 
frame; the field is canonically quantized on a three-dimensional space domain Q, 
with suitable boundary conditions on the frontier dfl. We are interested in the 
stress-energy tensor T^;^ or, more precisely, in the vacuum expectation value (VEV) 
of each stress-energy component: this is the so-called local problem in the theory of 
the Casimir effect. 

To deal with the divergences related to this problem, we introduce a "zeta regularized 
field" 0", depending on a complex parameter u and coinciding with for u = 0; 
this is used to define a "zeta regularized stress-energy tensor" T^^,, with finite VEV. 
The final step in this construction is the analytic continuation at m = of such 
VEV, which give the Casimir stress-energy for the case under consideration. A very 
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pleasant feaure of this approach, typical of the zeta regularization method, is that 
one gets directly a finite expression for the Casimir stress-energy, with no need to 
remove or subtract divergent terms. This is a major difference with respect to other 
renormalization schemes, employed more frequently for the Casimir effect; among 
these alternative approaches, let us mention the point splitting method, occasionally 
considered in this paper for a comparison. 

In the present work, the local zeta approach is mainly illustrated in the case of 
a field between two parallel plates, with Dirichlet boundary conditions (so, in this 
example fl = {—oo, +00)^ x (0, a), with a the distance between the plates); by simple 
variations of this setting, we then pass to the case of a field in the region outside 
one or two plates. The results obtained in these cases by local zeta regularization 
are compared with the ones derived by Milton [Sj and by Esposito et al [5] by point 
splitting, and they are found to coincide (incidentally, we take the occasion to show 
that the Casimir stress-energy tensors found in [8] and [5j, even though seemingly 
different, are in fact equal). 

To conclude this Introduction, let us briefly describe the organization of the paper. 
In Section [2] we sketch the basic framework for the Casimir effect, in the situa- 
tion outlined before: a canonically quantized scalar field on a space domain Q in 
Minkowski space-time, with given boundary conditions, its stress-energy tensor and 
the related VEV. In Section [3] we introduce the local zeta regularization scheme. 
In Section HJ we apply this scheme to the case of a Dirichlet field bewteen parallel 
plates, giving all details about the necessary analytic continuations; comparison is 
made with the results of [8] [5] . In Section [5l by simple geometric variations on the 
same theme, we derive the Casimir stress-energy tensor outside one plate, or two 
parallel plates. In Section |6l the outcomes of Sections HJ |5] are combined to derive 
the Casimir pressure on two parallel plates. 

For completeness, in Appendices \M and [B] we give some mathematical background 
on the analytic continuation of the polylogarithm, the function mainly involved in 
local zeta regularization for the case under investigation. 

2 Background for the scalar Casimir effect. 

Throughout this note we work in Minkowski space-time, which is identified with R"^ 
using a set of inertial coordinates 

X = (X'^)^=0,l,2,3 = {t, x\ x\ x'') = (t, X) . (2.1) 

We work in units where c = 1, = 1; the Minkowski metric is {rj^v) = diag(— 1, 1, 1, 1) 
(and is used to raise and lower indices). 

Let us fix a space domain Q C R'^ where we consider a neutral, massless scalar field 
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0; so, we have 

^■.Rxn^Csain) ; = 00 = (-a« + A)0 . (2.2) 

Here we are considering the space J0,{T-L) of hnear operators on the Fock space T-L, 
and the subset Cgai'H) of the selfadjoint operators; □ := d^^d^ is the d'Alembertian 
and A := Yl^=i 3-Laplacian. We assume appropriate boundary conditions 

(e.g., the Dirichlet conditions 0(t,x) = for x e dO). 

To expand the field in normal modes, we consider a complete orthonormal set 
{Fk)keK of eigenfunctions for the Laplacian in L^(Q,C), with the given bound- 
ary conditions; is a space of labels, for the moment unspecified, and we write the 
eigenvalues in the form —ujI- So, 

Fk-.n^C; AFfe = -ulFk (uk > 0); (2.3) 

/ d^^Tk{x)Fh(^) = 5(k, h) (k, heK) . 
Jn 

Any eigenvector label k can include different parameters, both discrete and con- 
tinuous. We generically write J^^dk to indicate summation over all labels, (i.e., 
literal summation over the discrete parameters and integration over the continuous 
parameters, with a suitable measure); 5{h.k) = 5{k.h) is the Dirac delta function 
for the labels space K (this reduces to the Kronecker symbol in the case of discrete 
parameters). The functions 

/fe : R X O ^ C , h{x) := Fk{^)e-'^>'' (2.4) 

fulfill n/fc = 0, and allow a unique expansion 

^^^'> = \akfk{x) + alfk{x)^ (2.5) 

(with ^ indicating the adjoint operator, and the complex conjugate). The destruc- 
tion and creation operators ak,al € ^{^) fulfill the relations 

[ofe, o/i] = , [ofc, ajj] = 5{h, k) , afc|0) = , (2.6) 
where |0) e is the vacuum state (of norm 1). 

Let us pass to the stress-energy tensor. This depends on a parameter ^ e R, and 
its components T^j, : R x Q — > Csai^H) are given by 

f^, := (1 - 20 d^^odj- Q - 2^^ rj^^d^^dj- 2^0 o 9^,0 ; (2.7) 
in the above, we use the symmetrized operator product Ao B := {1/2){AB + BA). 
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To be precise, a theory involving mereley a scalar field in flat space-time has a 
stress-energy tensor as above, with ^ = 0; the general form (12. 7p . with an arbitrary 
^, can be interpreted as the Minkowskian limit of the theory of a massless scalar 
field coupled with gravity, with ^ as the coupling constant, (p). 
Other authors have considered the general form fl2.7p independently of the previous 
interpretation in terms of a gravitational coupling; these authors invoke the principle 
that one can add to the stress-energy tensor a symmetric tensor with vanishing 
divergence, and regard the terms proportional to ^ in (12. 7p as additions of this kind 

m- 

To conclude these comments about ^, we mention that the choice ^ = 1/6 gives a 
conformally invariant theory [12], where the tensor (12. 7p has vanishing trace. For 
the above reasons, the term conformal coupling is usually employed to describe the 
case ^ = 1/6; one also speaks of a minimal coupling to indicate the case ^ = 0. 
After this digression, we proceed towards the Casimir effect considering the vacuum 
expectation value (VEV) of T^^^. We use the expansion (12. 5 p for the field, with the 
relations {0\akah\0) = 0, {0\alal\0) = 0, {0\alah\0) = and (0|afca[|0) = 6{k,h)] in 
this way, we readily obtain the formal expression 

(i ~ 2 ) i^^f'^'^'f'' + ^'^/'^^M^) (2-8) 

however, the above integral is divergent and some renormalization procedure is 
needed. A standard approach relies on the so-called point splitting regularization 
(see [21 131 [5]; [8] essentially uses the same method). In this approach, in place of 
T^y{x) one considers 

^In the theory of a classical, massless scalar field coupled with gravity, the dynamical vari- 
ables are the field (j) and the space-time metric g^^. The action functional is S[(j),g^i,] = 

i / d'^x^—g(j)'^(f)df_i(t> — i?(g^ — in units where the gravitational constant is 1, where 

g := det{gf_i^) and R is the scalar curvature of the metric: for more details see, e.g., [12] page 
43. One imposes the stationarity of the action with respect to variations of 4> and gf_iu- in this way 
one gets, respectively, the scalar field equation and Einstein's equations with the field stress-energy 
tensor T^^. One can analyze the almost Minkowskian case where (p is small and 5^1/ = 77^1/ -I- /i^^iy, 
with h^i, a small perturbation of the second order in (j). In this case the scalar field equation has 
the form D(f> — 02[<p] with □ the Minkowski d'Alembertian d^d^^ and the stress-energy tensor is 
T^^v = (1 - 2^) d^(j)d^(j) - - 2^) r]f^^d^(j)dx(t> - '2-£.4>dtii,(l) + 02,[(j)\\ here, O2 and O3 indicate terms 
of orders 2 and 3. Of course, Einstein's equations relate h^^ to (p. Neglecting the higher order 
terms, and quantizing the field, we obtain Eq.s (12. 2p (|2.7p . 



(o|r,.|o) = / ^ 

J K 
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%,{x,x') ■= (1 - 20 d,^{x) o - Q - 2^) v,ud^^{x) o (2.9) 

-2^ 0(x) o (9/,^0(x') , (x, a;' G R X (]) , 
giving formally T^^[x) in the limit x' — )■ x. One then defines the renormalized VEV 

(0|f^,(x)|0),e„ := FP (0|7;,(x,x')|0) , (2.10) 



of T^y{x) as 



where we have written FP to indicate the "finite part" in the limit x' — )■ x; this 
means that one writes down the VEV of T^,y(x, x') and then removes the terms 



diverging for x' — )• x. (jj) 
Hereafter we will describe the alternative approach considered in this paper, i.e., the 
local zeta method. 



3 Local zeta regularization. 

In the sequel we keep all the notations of the previous section. Let us denote with u a 
complex parameter and consider the powers (— A)~"/'^, built from the 3-dimensional 
Laplacian A. From AF^, = —u\Fk it follows (— A)~"/^Ffc = w^^^^F^, whence 

(-A)^"/V. = c^r^V. ; (3.1) 

there are similar relations for the conjugate functions, starting from AF^ = —ulFk- 
We now introduce the smeared, or zeta-regularized field operators and stress-energy 
tensor 

0" := (-A)-"/^0 , (3.2) 
f;, := (1 - 20 d,^^ o dj- - Q - 2^^ v,.d^rdxr - 2^ 0" o 8,,^ , (3.3) 
which formally give and T^i, in the limit m — i- 0. Eq. (12. 5p implies 

akfk{x) + a\fk{x) ; (3.4) 



X 



/K u/2+1/2 



^Oi course, the concept of "finite part" contains a basic ambiguity, that must be removed by a 
precise prescription. In the case of an electromagnetic field in Minkowski space-time, a prescription 
of this type has been given in [2]; this approach could be adapted to the scalar case. An alternative 
strategy is to define the finite part in (|2.10p as the x' ^ x limit of what remains after substracting 
from (0|T^ii/(a;, a;')|0) the analogous VEV for a field without boundary conditions. For a critical 
analysis about these and other problematic aspects of point splitting, see [T^ . 
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now, a computation very similar to the one giving Eq. fl2.8p produces the result 



(0|T"|0) 



\~\)[ 9f^fkdufk + dufkdf^fk 



~{\~ ^^i^^^^fkdxfk - ^ ( fkd^^ufk + fkd^^ufi 



The above integral typically converges for 3ftu sufficiently large and is an analytic 
function of u, a situation that will be exemplified hereafter. Eq. (13. 5p . with 3ftu 
sufficiently large, is our regularization of the VEV for T^^] we now define the renor- 
malized VEV as 



(0|T^.|0), 



AC 



u=0 



(0|T"|0) 



(3.6) 



indicates that one should consider the analytic continuation of the 
|T" |0), and evaluate it at m = 0. In the next section, the whole 



where AC 

u=0 

function u i— )■ 

procedure will be exemplified in the classical case where Q is the region between 
two parallel plates, with Dirichlet boundary conditions; in the subsequent section 
we will treat the region outside one or two plates. 



4 Casimir effect between two parallel plates 

Setting up the problem; the zeta-regularized stress-energy tensor. Let the 

plates occupy the planes = and = a (a > 0); the region between the plates 
is 

n := {{x\x^, x^) \ x\x^ eR ,0 < x^ < a} . (4.1) 
We assume the Dirichlet boundary conditions 

0(^,x^x^x^) = forx3 = 0,a. (4.2) 

Let us produce a complete orthonormal set {Fk)k&K of Dirichlet eigenf unctions for 
A on fi, and the corresponding eigenvalues We can take 

K := {k = {ki, k2, ks) \ ki, k2 ^ H , k^ E {Tv/a, 27r/a, Sir/a, ...}} , (4.3) 



1 



dk := / dki / dk2 



A:3€ {7r/a,27r/a,37r/a,...} 



7rv2a 



The above functions fulffil d^^FkFh = S{ki — hi)5{k2 — /i2)^fc3,/i3; we can use them 
to build fk(t,:s.) := Ffc(x)e~*'^''*. Let us pass to the computation of the components 
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(0|T^^|0), and to their analytic continuation at m = 0; we will start from the case 



/i = 0, z/ = 0. From Eq.s ([33]) and fl4.3tl4.4l) . we obtain 

(0|To"o|0) 



(4.5) 



1 



Svr^a 



E 



fcsG {7r/a,27r/a,37r/a,...} 



R 



dki / dh 



^ +00 „ „ 



kj + e^ + kj- {kj + A;| + A^kj) cos{2ksx^ 
R ' {kl + kl + fci)V2+i/2 

g2 + g2 + £2 _ (^2 ^ ^2 ^ 4^£2) cos(27r£xVa) 
' {ql + g2 + ^2)«/2+i/2 



where, in the last passage, we have performed a change of variables ki = (7r/a)gi, 
^2 = (7r/a)g25 ^3 = {T!'/a)i. We now pass to polar coordinates in the {qi,q2) plane, 
setting qi = pcos9, q2 = psm9, and then compute the integrals in 9,p; in this way 
we obtain 



oo /■2tt 



(0|To"o|0) 



£=1 ''^ 



dO 



+ 00 



dpp 



1 



47r"-2(M - 3)a4 



£=1 



1 



p^ + e- (p^ + 4^^^) cos(27r£xVa) 

(p2 _^ £2)«/2+l/2 

2 + 4(m - 3)^ cos(27r£xVa) 



{u-D 



£n-3 



The last series is clearly convergent if 



5Ru > 4 ; 



• (4.6) 



(4.7) 



under the same condition, all the expressions given previously for (0|Tf5{)|0) are mean- 
ingful and finite. To go on let us recall that the polylogarithm (-2, s) i— )■ Lis{z) is 
defined by 



+ 00 / 



Lis{z) := — for z G C, \z\ ^ 1 and s G Sz 



where 5^ C C is the set of values of s for which the above series converges: one finds 

(4.9) 



5. 



C ifkl<l, 

{^s > 0} if 1^1 = l,z=^l, 



{3fts>l} if 2 = 1. 



(We note that, for \z\ > 1, there is no s G C such that the series converges.) Let us 
also recall that the Riemann zeta function s ^-^ (^{s) is defined setting 



+ 00 

C(5):=Lz,(l) = 5^- for.GC, 3?.>1 . 



(4.10) 
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The functions Li, ^ can be extended to larger domains by analytic continuation, 
as reviewed in Appendix |Al Comparing Eq.s (14. 6 p (14.81) (I4.10p . and noting that 
cos(27r£xVa) = (1/2) (e^^^^'/'*)^ + (1/2) (e-^^^^'/'^)^ we see that 



1 



47r"-2(M - 3)a4-« 



C(«-3)- 



(0|To"o|0) 
l + 2(M-3)e 



Liu-3{e 



+ Liu-sie 



(4.11) 



for > 4. The other components (0|T^^|0) are treated similarly. More precisely, 
we find 

mm 



1 



E 

ks&{iT/a,2iT/a,3n/a,...} 



R 



dki / dh 



R 



A;2 - {kf + (1 - 40A;i) cos(2A;3x3 
' {kf + kl + A;|)«/2+i/2 



(01^3-310) 



(4.12) 
= 1,2); 

(4.13) 



fcsG {7r/a,27r/a,37r/a,...} 



R 



(i/ci / dk 



K ^(A;2 + fc2 + A;|)«/2+i/2 



(o|r"|o) = o ^oI^^^v 



(4.14) 



indeed, one checks that (0|T^2|0) = with a change of variables k2 ^ ki. The 

expressions (I4.12p (I4.13P can now be treated with the same method employed for 
(0|T(]^|0): one makes a change of variables ki = {n/a)qi, k2 = (vr/a)g25 ^3 = (^r/ct)^, 
passes to polar coordinates (p, 9) in the (gi, ^2) plane, integrates in these coordinates 
and then expresses the remaining sum over i in terms of the zeta function and of the 
polylogarithm. The results of such computations can be summarized in the formula 



r"|o) 



At,i^=0,l,2,3 



+5"(a;3) 



/ -i-2(M-3)e 

1 




C(^-3) 
47r"-2(M-3)(M- l)a4-«^ 



u 
2 



/ M- 1 \ 

10 

10 

\ u-3 ) 



'2(n-3)^ 





(4.15) 



0\ 









l-- + 2(M-3)e 




for "Siu > 4, 



0/ 



-2inx^ /a 



47r"-2(M-3)(M-l)a 



A-u 
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Renormalization by analytic continuation. Due to fl4.15p . the problem of the 
analytic continuation of (0|T^^|0) at m = is reduced to the problem of continuing 
the functions s t-)- ({s), Lis{z) (for fixed z) up to the point s = —3. As reviewed in 
Appendix [B], such continuations are given by 



z{z^ + 4^ + 1) 



for;.^!; (M) = L^„3(l) = ^ 



(4.16) 



(Note a discontinuity with respect to z presented by the continuations: lim^^i Li^3{z) 
= oo 7^ Li_3(l). For an interpretation of this fact, we refer again to Appendix [B|) . 
Keeping in mind these facts we return to Eq. fl4.15p . from which we infer that 



(0|T^.|0), 



(0|T^.|0), 



AC 



u=0 



T" |0) is as follows: 



= A 



-10 

1 

1 

\ 



\ 




-3/ 



/ -1 \ 

10 

10 

\ / 



(4.17) 



TT 



HereA:=A°: 

12a4 

i.e., using the expressions (14.16^ 



<(-3) andS(x3):=B0(x3): 



TT 



Li-3{e 



2inx^ /ay 



'2iTTX^ /a 



A 



1440a4 



(4.18) 



B(x^ 



2 + cos(27rx^/a) vr^ 3 - 2 sin^(7rx^/a) 



12a4 [1 - cos(27rxVa)]2 48a4 sin^(7rxVa) 



(0 < a;^ < a) 



(The first expression above for B[x^) follows using ( Km with z = e±2*^^'/"; the 
second expression follows from the duplication formula for the cosine). 
Let us remark the following: 

i) For /i = 0,1,2 the components (0|T^^|0)ren depend on x^ through the function 
B, except in the conformal case = 1/6 where they are constant. The component 
(0|T33|0)ren IS coustaut iu any case. 

ii) The function B{x^) diverges like l/(x^)^ in the limit x^ — )■ 0, and like l/{x^ — a)^ 
in the limit x^ — )■ a. The same can be said of (0|T^^|0)ren for /i = 0, 1, 2 and ,^ 7^ 1/6. 
Eq.s f l4.17lH?T5]) are our final result for the renormalized stress-energy VEV. We have 
now checked the following claim of the Introduction: the local zeta method, based 
on analytic continuation, gives directly a finite stress-energy tensor, with no need 
to remove divergent terms. We already indicated this fact as a relevant difference 
between this approach and the point splitting method; however the renormalized 
tensors derived by these two approaches coincide, as illustrated hereafter. 
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Comparison with the results obtained by point sphtting. Let us compare 
our Eq.s fl4.17ti4.18|) with the resuhs obtained by Esposito et al. [5J by the point 
splitting method. The essence of this method has been reviewed in Eq.s (12.9112. 10]) 
(which are implemented in [5] using a Green function method, fully equivalent to 
the eigenfunction expansion for the Laplacian). The cited work produces the formal 
result 



lim (0|f^,(a;,x')|0) 



A + 



lim 



27r2 ^3'^^3 (a;3 - x^'Y 

/ -1 \ 





-1 








\ 







1 
















1 





V 











-3^ 



(4.19) 



10 
10 
\ / 



where A, B are as in our Eq.s f l4.18p . As indicated in fl2.10p . in this apprach renormal- 
ization amounts to subtract the divergent term proportional to Ymv^-ii ^^'i{x^ — x^')~^] 
so, the renormalized stress-energy VEV agrees with ours. 

A stress-energy VEV renormalizazion, based essentially on point splitting, appears 
as well in the previous book of Milton [8] who gives for A the expression in (I4.18P 
but obtains, in place of S, the function 



'Bix-' 



[C(4,xVa) + C(4,l-xVa)] 



167r2a4 

here (z, s) (-)■ C(s, z) is the Hurwitz zeta function defined by 

+00 



(4.20) 



c(^,^) = E 



1 



1=0 



(4.21) 



Indeed, the Milton function S coincides with the function B in Eq. f l4.18p . To show 
this, we refer to the known identity (see [H], page 608, Eq. (25.11.12)) 



as+i,z) 



■1) 



ij^'\z) for s = 1,2,3,... , 



s\ 



(4.22) 



where the right hand side contains the polygamma function ip^'^^ (z) := (d/dz)'^^^ In r(z) 
for s = 1, 2, 3, this implies 



S(a;^) 



1 



[ij^''\xya) + ^P^^\l-xya)] 



(4.23) 
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Another relation, known to hold for the polygamma function, is 



^P^'\l-z) + {-iy^^^^'\z) = {-iyTT^cot{Trz) for s = 1,2,3,... ; (4.24) 



(see [TT], page 144, Eq. (5.15.6)); this entails 



COt{lTz) 



TT^ [3 - 2sin2(7ra;Va) 



(4.25) 



48a4 



sin^(7rx3/a) 



B{r') as in f gl8|l . 



As a final comment on this result, we mention that the equality !B = i? is a special 
case of a more general relation between the polylogarithm and the Hurwitz zeta 
function (see f7], or [TTj for a reformulation in modern notations). 



5 The Casimir effect outside one plate, or two 
parallel plates. 

The case of a single plate. Let the plate occupy the plane = 0; heferafter we 
determine the renormalized VEV of T^jy in one of the half-spaces bounded by the 
plane, say, in 

fioo := {{x\ x^, x^) I x\ e R , > 0} . (5.1) 

As before, we assume for the (scalar, gravity coupled) field (p the Dirichlet boundary 
conditions 

${t, x\ x^, x^) = for = . (5.2) 

To treat this case, it is not even necessary to set up a framework as in the previous 
sections, starting from the Dirichlet eigenf unctions of A in Q^o- In fact, it suffices 
to view Qoo as the a — )■ +oo limit of the domain 

Qa ■■= {{x\ x^, x^) \ x\x^ eR ,0 < x^ < a} (5.3) 

and define the renormalized VEV of T^,y in f2oo as 

oo, ren ■ 

lim (0|T^,|0) 

a, ren j 

(5.4) 

a— >'+oo 

where the right hand side contains the renormalized VEV in Qa', the latter is known 
from the previous section, see Eq.s (14. 17114. 18|1 . So, 
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\T,u\0), 



lim Ac 

a— >+oo 



/ -1 \ 

10 

10 

\ -3 / 



(5.5) 



+ (1-60 lim Ba{x^ 



/ -1 \ 

10 

10 

\ / 



, Aa := A, Ba{x^) := 5(x3) as in dHH]). 



From Eq. ( Hl8|) . it is evident that A„ -> 0, ^^(x^) ^ l/(167r2(x3)^) for a +oo; 
so, 



(o|r^.|o)oo,. 



1-6^ 
167r2(x3)4 



/ -1 \ 

10 

10 

\ / 



(0 < < +oo) . (5.6) 



Of course, one 



The above result, derived in a different way, appears e.g. in 
obtains similar conclusions in the half space {— oo < < 0}. 
We observe that, if the coupling parameter takes the conformal value ^ 
stress-energy tensor vanishes everywhere outside the plate. 



1/6, the 



The case outside two parallel plates. We now consider, as in the previous 



section, two plates occupying the planes x = and x 



we are interested in 



the renormalized VEV of T^j, in the region outside the plates, which is the disjoint 
union of the half spaces {x^ < 0} and {x'^ > a}. This can be obtained by obvious 
adaptations of the result (15. 6p on the half space {x^ > 0}; the conclusion is 



(o|r^.|o), 



(0|T^.|0),„ 



167r2(x3)4 



1-6^ 



167r2(x3 - a)4 



/ -1 \ 

10 

10 

\ / 



/ -1 \ 

10 

10 

\ / 



-oo < < 0) . 



(5.7) 



(a < X < +oo) 



(5. 



Note that, in the conformal case ^ = 1/6, (0|T^^|0)ren is identically zero outside the 
plates. If 7^ 1/6, the components with /i = z/ = 0, 1, 2 of this tensor diverge like 
l/(x^)^ and l/(x'^ — a)^ for x'^ — > 0^ and x^ — )■ a"*", respectively; we recall that similar 
divergences were found as well for the stress-energy tensor between the plates. 
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6 Pressure on the plates 



In this section we alway use the spatial indices i,j e {1, 2, 3}. Let us consider any 
one of the two plates ai = or = a, and evaluate the force per unit area 
acting on it; in principle, this computation should take into account the action of 
the field both inside and outside the plates. The force per unit area produced on the 
given plate by the field in the inner region is = (0|T*j |0)m^o„( where (0|T*^ |0)m 
is the renormalized stress-energy tensor in the inner region and n^^^ the normal 
unit vector to the plate pointing towards the outer region. On the other hand, the 
force per unit area produced on the same plate by the field in the outer region is 
Pout — (Ol^MO)out^^iri' where the subscripts j„, out have an obvious meaning. So, the 



^out \"|-^ jl^/ouf'in^ 

total force per unit area on the plate is 



P 



in|0).nnL+(0|r,|0),„,n^'„. 



(6.1) 



For the plate located at x^ = 0, we have (nl^) = (0, 0, 1 



in- 



out) 



(0,0,-1), so 



p' 



(o|mo)„„,-(o|r3|o), 



for the plate at x^ = a the inner ond outer normals are reverted, so 



(6.2) 



-{0\T\\0)out+{0\T\\0)i 



(6.3) 



Now, we take the expressions of {0\TWO)in,out = (0|Tj3|0)i„_out from Eq.s (14. 17114. 18|) 
and (OESD; for both plates (0|fi3|0)o„t vanishes and ((0|fi3|0)i„) = (0,0, -3A) 
with A = 7r^/1440a^, as usually; in conclusion 



P 



(0,0, 



(p' 3 ) =(0,0, 



480a4' 
7r2 



(6.4) 



(6.5) 



480a4' 

Thus the plates are subject to a reciprocal attraction inverserly proportional to the 
fourth power of their distance. We note that, once more, the result obtained agrees 
with the ones reported in [5[ |8]. 
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A Appendix. Analytic continuation of the poly- 
logarithm (and of the zeta function). 



Let us report Eq.s (gSD fICTD 

Lis{z) := — ior z E C, \z\ ^ 1 and s G Sz] 



1=1 



+00 

Cis) := Ltsil) := 5^ - for s G C, 3fts > 1 



E=l 



(In the above 5*2 is the subset of C such that the series for Lis (z) converges, see Eq. 
(14. 9p ). Our problem is continuing analytically (in s) the functions defined as above; 
the solution is well known [7], and reported here for completeness. Indeed, let us 
define 



T(l 



s] z 



2m 



dt 



'-t) 



s-l 



e- — z 



for ;s G C, s G C \ {1,2,3,..} , (A.l) 



Lis{z) := lim Lv(z) for 2 G C \ {1}, s G {1, 2, 3, ...} or z = 1, s G {2, 3, ...}; (A.2) 

C(s):=Li,(l) forsGC\{l}. (A.3) 
In Eq. (lA.ip . F is the usual Gamma function. Furthermore: 



i) Hz is a Hankel contour in the complex t plane, starting from infinity in the direc- 
tion of the positive real axis, turning counterclockwise around t = and returning to 
infinity in the direction of the positive real axix (see the figure below); this contour 
is chosen so that all the solutions t of the equation e* = z are outside the region 
bounded by Hz-, except the solution t = appearing if 2; = 1. 

ii) For each t E Hz v/e intend 



-t) 



s~l 



-i(s— l)7r^s— 1 



t 



s-l 



U|s-lgi{s-l) argt 



(A.4) 



where t i— )■ argt is the unique continuous function on Hz such that argt — )■ when 
t tends to the beginning of the path. 
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For each fixed z with \z\ ^ 1, the function s i— Lis{z) defined via flA.ip f lA.2p is the 
(unique) analytic continuation of the function s E Sz ^ Lis{z) previously defined 
via the power series (14. 8p : to prove this, it suffices to prove that the definition (lA.ip 



for Lis{z) via a contour integral implies a series expansion as in (14. Sp . if s G Sz- To 
this purpose, we reexpress the function of z and t in (lA.ip in the following way: 



z 



— z 1 — ze 



'-ty-'-r^^, = i-ty-' J2^^e-r (a.5) 



=1 



inserting this result into Eq. ( lA.ip . we obtain 



1=1 ^ 

On the other hand, the known Hankel's integral representation for l/F [11] implies 

for I = 1,2,3... . The last two equations yield the wanted expansion Lis{z) = 
J2t=i ^7^', of the form fOD . 

The above manipulations have hidden a problem: to grant convergence of the series 
expansion ( lA.Sp and the exchange between the summation over i and the integration 
over Hz, one shoud have \ze~^\ < 1 uniformly in t G -ff^: on the other hand, \ze~^\ 
can be larger than 1 when t is on the arc in the half plane 3fJt < 0, turning around 
the origin. Let us skecth how to overcome this difficulty; the basic idea is that 
Lis{z) defined in (lA.ip does not change if we shrink the path Hz around the origin. 



If \z\ < 1, we can shrink Hz so that \ze~^\ < 1 uniformly in Hz, including the arc 
that turns around the origin. The case = 1 is a bit more technical: one isolates 
from the integral over Hz the contribution of the arc encircling the origin, makes 
a series expansion of the integrand in the remaining part of Hz, and finally proves 
that the contribution from the arc can be made arbitrarily small by shrinking. 
The previous results on the analytic continuation of the function s i— )■ Lis{z) hold, 
in particular, for z = 1; so, the function s G C \ {1} i— )■ ({s) in flA.Sp is the (unique) 
analytic continuation of the function defined previously by (14.101) . 
Another property of the function (lA.ip(]A.2p is that it is jointly analytic in {z,s), 
when these variables range in a suitable open subset of C^; outside this open set, 
some pathologies can appear. In particular, for a given s, this function can happen 
to be discontinuous in z at the specific point z = 1: see, e.g., the case s = —3 
discussed hereafter. 
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B Appendix. The polylogarithm (and the zeta 
function) at s = —3. 

Let us consider the analytic continuation of the polylogarithm described in Appendix 
IX| and evaluate it at s = —3. For this choice of s, Eq.s (lA.ip takes the form 



Li. 



6 z 



dt 



2m t^{e' - z) 



for 2 G C : 



(B.l) 



the integral therein is easily computed by the method of residues, as briefly skecthed 
hereafter. First of all, the integral in ( lA.ip involves a meromorphic function of t, 
whose only singularity in the region bounded by is a pole at t = 0. The order of 
the pole is 4 if 2 7^ 1, while it is 5 if 2; = 1, and one finds 



Res 



[z^+Az + l] 



if z 7^ 1; Res 



i=0 



6(^-1)4 

When these results are inserted into fIB.ip . the residue theorem gives 

Ziz"^ +4:Z+l 



t=o 



720 



Li_^{z) 



[z-l] 



for^eC\{l}; a-^) = Li-?X^) 



1 
120 



these are the statements (I4.16p . which are now justified. 

To conclude, we note a discontinuity of the type mentioned at the end of Appendix 
Rl lim2_j.i Lz_3(z) = 00 7^ Li_3(l). This is basically due to the jump in the order of 
the pole (from 4 to 5) when z goes to 1. 
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